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Abstract: We characterize the two-site marginals of exchangeable states of a system 
of quantum spins in terms of a simple positivity condition. This result is used in two 
applications. We first show that the distance between two-site marginals of permutation 
invariant states on N spins and exchangeable states is of order 1/N. The second application 
relates the mean ground state energy of a mean-field model of composite spins interacting 
through a product pair interaction with the mean ground state energies of the components. 



1 Introduction 

The mean-field approximation is a very common approach in statistical mechanics. It 
consists in replacing suitably chosen parts of the interaction by their expectation values. 
This generally simplifies the problem but leads to non-linear self-consistent equations for 
the dynamics and the equilibrium states. This kind of approximation often leads to rea- 
sonable results in regimes where the interactions are rather weak. In other cases, the 
self-consistency equations may induce artificial phase transitions [11]. 

A characteristic feature of the most basic version of the approximation is that every 
particle interacts in the same way with every other particle. Therefore the Hamiltonian 
and the ground and equilibrium states have a huge symmetry: particles can be arbitrarily 
permuted. By mean-field models we here mean quantum spin systems which exhibit this 
kind of symmetry. There is a vast literature on the subject dealing both with the structure 
of states and dynamical maps [3, 7, 12]. We briefly recall some essential notions and results. 

A state uj of a system of iV identical spin-(d — 1)/2 particles is determined by a den- 
sity matrix p £ A4 d N(C) = <2>N-Md(C), where Afd(C) denotes the complex matrices of 
dimension d 

uj(A)=TrpA, for A £ ® N M d (C). 

An iV-particle state is symmetric if it is invariant under permutations of the particles, i.e., 
if 

f4(|*i)®---®|*jv)) = |* w (i))®---®|* w (j\r)> where |^)£C d , (1) 

then 

lo(A) = ^(U^AU*) for every permutation ir of {1, . . . , N} and A £ <8>jv.M<f(C). 
In terms of the density matrix p of u, 

P = u; P u 7r . 
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Consider a system of (N + M) particles with a symmetric state uj. For each subsystem 
of N particles the marginals of uj are symmetric iV-particle states uj^ 

uj n (A) :=uj(A® (® M 1)) for every A E ® N M d (C). 

Note that, due to the symmetry of uj, only the number of spins in a subsystem matters 
and not the precise sites on which the subsystem lives. The density matrices pn associated 
with these states are obtained by taking partial traces of the density matrix p that defines 

UJ 

p N = Tr M p= ( id ®l e u •••e iM )(e il ■■■e iM |) (p), 

(u,...,im) 

where {ej}f~Q is a basis of C d . The converse is not true, a symmetric iV-particle state uj 
cannot always be extended to a symmetric (N + M)-particle state. For example, consider 
the pure two-qubit state determined by | ^) = 4=(| 01) + | 10)). This state is symmetric 
but has no symmetric extension to three qubits [13]. 

If we want a symmetric state to have a symmetric extension to an arbitrarily large sys- 
tem, we have to impose the stronger condition of exchangeability. A state uj on ®jv-M<i(C) 
is called exchangeable if it admits for any M > a symmetric extension uj(n+m) to 
®n+m •Md(C). Exchangeability is a quite strong condition, as we see in the following 
quantum version of de Finetti's theorem [1, 6]. 

Theorem 1. If uj is an exchangeable state on <S>jv-Md(C), then 

uj — dp(a) <g)jv a 
Js d 

where Sd denotes the state space o/A4d(C) and p is a probability measure on Sd- 

The exchangeable states are mixtures of symmetric product states which implies that 
they are non-entangled and so only classical correlations are possible [10]. The inverse 
implication is not true, not every symmetric separable state is exchangeable. Consider, for 
instance, two density matrices p, a e j\4d(C), then the state associated with |(p®o" + (T®p) 
is symmetric and separable on A4d(C) <E> j\4d(C) but generally not exchangeable. 

2 Two-site marginals of exchangeable states 

We want to characterize the exchangeable states on two particle systems with d degrees of 
freedom. 

Theorem 2. A symmetric state uj on .M^C) <8> M.d{C) is exchangeable iff 

uj(B ®B)>0 for all B G M l d (C) 
where A4^(C) denotes the complex Hermitian matrices of dimension d. 
Proof. If uj is an exchangeable state two-particle state, then by theorem 1 we have that 

u{B ®B)= I dp{a) a{Bf > for every B E M d (C). 
Js d 

The remaining of the proof is postponed until section 2.2. □ 



2 



In order to prove the inverse direction we use the polar cone theorem to invert the role 
of states and observables. So, instead of proving that u is exchangeable if uo{B <g> B) > 
for every B G Ai^(C), we will prove that a flip-invariant, hermitian A G -M<f(C) <8> A^^(C) 
is a positive combination of B a <g> £? Q , _B a G .Mj^C), if Tr(<r <g> a A) > for every density 
matrix a G Sa- 

More explicitly, given a real Hilbert space H and a set C C 7Y, the cone 

C* := | (x,y) > for every x G C}, 
is called the polar cone of C. 

Theorem 3. Let H be a real Hilbert space and C a subset ofH, then 

(C*)* = Cone(C), 
where Cone(C) denotes the closure of the cone generated by C . 
Let F be the flip operator on C d <g> C d 

F(<p ® ip) := ip ® (p. 

We consider the real subspace /C of the complex hermitian matrices of dimension d 2 which 
commute with F and equip /C with the trace scalar product 

(•, •) : M h d2 (C) x M)p{C) -> C : (At, A 2 ) ^ Tr A 1 A 2 . 

We choose C to be the set of all symmetric two-site product states determined by density 
matrices on C d 

C = {p (g> p | p G A^rf(C), p is a density matrix}. 
It is then enough to prove that the polar cone of C is the closed cone C* generated by 

{B®B | BeM h d {C)} U {L | Le (M d (C)®M d (C)f, L > and LF = F Lj 

where L > means that L is a positive semi-definite matrix. Indeed, applying the polar 
cone theorem, we get 

C** = {p | Trp(B ® 5) > 0, B G 7Wrf(C) and TrpL > 0, L > 0} 
= Cone(C) = Cone({p ® p | p G «S d })- 

We shall first prove the analogous result for the classical case, that is when we replace 
the matrix algebra A^rf(C) by the diagonal matrices of dimension d and states by probability 
measures on the relevant configuration space. 

2.1 A classical intermezzo 

Let Q be a finite set and 

C := {p x p \ p is a probability measure on Q} 

then 

C* = {/ : ft x ft - R = and 

(p x p)(f) > for all measures p on ft}. 

The aim is to show that the cone C* is generated by functions of the form f\ + f 2 where 
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i) fa > and fi(x,y) = f±(y,x) 

ii) f2 = 9 x 9 with g : Q — > R 

Fix / in the interior of C*. By subtracting from / a suitably chosen non- negative 
symmetric function, we can arrange to have a strictly positive measure p on Q such that 

(/i x /io)(/) = and (/i x //)(/) > for all measures //. (2) 

Let p now be a measure as in (2). For any r, a sufficiently small real functional on Q, 
p + t is non-negative on f2. Therefore, by assumption, 

((p + r )x(po + r))(/)>0. (3) 
As (p x po)(/) = 0, this can only hold if 

(/i x r)(/) = for all choices of r on H 
In this case, condition (3) translates into 

(rxr)(/)>0, for all r. (4) 

As the matrix F := [f(x,y)] is real and equal to its transpose, (4) amounts to requiring 
that F be semi-definite positive. But then there exist Cj(x) such that 

f(x,y) = [F] X:V = ^[^-(a:)^-^)], 
j 

proving our statement. 

2.2 The quantum case 

Proof of second part of Theorem 2. We have now C := {p <E> p \ p is a density matrix in 
M d (C)} and 

C* := {A e M%(C) \ AF = F A and Tr A(p® p) > V density matrices p G M d (C)}. 

The aim is to prove that the cone C* is generated by matrices of the form A 1 + A 2 with 

i) Ai > and A 1 F = F A x . 

ii) A 2 = B®B with 5 G M$(C). 

As in the previous section we fix A in the interior of C* and subtract from A a positive 
semi-definite matrix to have an invertible density matrix p such that 

Tr A(p <S> Po) = and Tr A(p <g> p) > for all density matrices p. 

For any choice of B G A4^(C), with sufficiently small, po + -B is still positive semi- 
definite and so 

Ti A((p + B) ® (p + B)) >0. (5) 
As Tr A(p ® po) = 0, this can only hold if 

Tr(p ® fi)A = for every 5 G A<5(C). 
4 



In this case, condition (5) translates into 

Tr(B® B)A > for every B e M d (C). (6) 

We now extend the argument for the classical case, see section 2.1 to the quantum case. 
Therefore we introduce real linear maps, 

V d :M^(C)^H and 

M d : (M d (C)®M d (C)) h ^B(H), 

where H is a suitably chosen real Hilbert space and B{7i) denotes the linear operators on 
that space such that: 

i) V d and M d are one-to-one and onto. 

ii) For every A G C*, M d (A) is positive, this will follow from condition(6), M d (A) T = 
M d (A) and Tr A(B ® B) = (V d (B) \M d (A)\ V d {B)). 

iii) M?(\t)(t\)=V;\t)®V;\t). 

With these maps we can prove that A = B a ® B a . Indeed, as 

Tr A(B <g> B) = (V d (B) \M d {A)\ V d {B)) > for every B E M h d (C) 

and V d is onto, we get that M d (A) is positive or M d {A) = J2 a \ T a ){r a |. Now, because M d 
is one-to-one and using property (iii) above, we have 

A = M ^ (E I ^ l) = E M d 1 (\ T«)(Ta |) = E ^ _ Va) ® ^K), 
ol a a 

proving our statement. Constructing the maps V d and M d and verifying their properties is 
rather tedious. We therefore provide the details separately in appendices A-C. □ 



3 Finite size symmetric states 

In this section we focus on the distance between the two-site marginal of an iV-particle 
symmetric state and the two-site exchangeable states. Let S N be the set of symmetric 
states u) of two particles which have a symmetric extensions to iV sites and let S°° be the 
exchangeable two-particle states. Obviously, 

S 2 dS 3 ---dS n dS n+1 ---dS°°. 

The sets S N are closed and convex in the state space of Ai d (C)® M. d (C) for all N — 2, 3, 

We can now wonder about the distance of S N to the exchangeable states <S°° 

d(S N ,S°°) = maxd(w,5°°) = max min ||w-w'|| 

= max min Trip — p'\, (7) 

where p and p' are the density matrices corresponding to the two-site states oo and u'. 
We know that for N — > oo, this distance vanishes, but we are interested in the behaviour 
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with N. An upper bound of the order 1/vA was obtained in [9]. Such bounds yield a 
measure of the maximal entanglement of states in S N . For a detailed analysis of a model, 
see e.g. [2]. 

A possible approach to this question is to use the information on the structure of 
symmetric states that can be obtained from group theory. The decomposition of the natural 
representation of the permutation group Sn of a set of A elements on (C d )® N given in (1) 
in irreducible representations is a highly non-trivial achievement of group theory [5]. The 
result is that the irreducible representation of Sjy are labeled by standard Young tableaux 
T. The irreducible representation corresponding to T has dimension d(T) and occurs with 
a multiplicity m(T), both d and m are explicitly known, moreover, d depends on A and 
m on A and d. Hence, there is a decomposition 

(C i )* N = eC B ^®C^. (8) 

T 

Any symmetric A-particle density matrix is then of the form 

p = c(T) p T <g> 1 (9) 

T 

where pr is a density matrix on C m( - T ) and c(T) are suitably chosen non-negative normal- 
ization coefficients. In order to estimate the distance (7) we can compute the two-site 
marginals of a state determined by a pure in (9) and estimate its distance from the 
exchangeable states. Such a computation is, however, rather involved. We nevertheless 
sketch an example of the computation for the case d — 2. 

Considering C 2 as the state space of a single spin-1/2 particle, the decomposition (8) is 
nothing else than the standard decomposition of a system of A spin-1/2 particles according 
to total spin. Any value of the spin in {0, 1, ... , A/2} for even A and {1/2, 3/2, . . . , A/2} 
for odd A occurs. Let us simplify the problem even further by choosing a completely 
symmetric normalized vector \1> in (C 2 )® N . We fix canonical basis vectors | |) and | j) in 
C 2 , e.g. to the eigenstates of the z-component of the spin. A natural basis of the completely 
symmetric subspace of (C 2 ) 81 is then { \n) | n — 0, 1, . . . A} where \n) is the normalized 
state obtained by symmetrizing an elementary tensor with n factors | f ) and A — n factors 
||). Our vector ^ can then be written as 

N 



V = J2 a n\n), (10) 



n=0 

where the a n are components of a normalized vector in C N+1 . To calculate (*ff \X\ ^), we 
need to know the (m \X \ n). We are especially interested in 

X = A e M 2 and X = M e M 2 ® M 2 - 

A possible trick is to consider 

TV 2 N 

x = ® N e sA = i + j2^ + s -( Yl A i® A i+ J E A *) + --- 

i =1 {i,j\i¥=j} i =1 

with A G M 2 - Then 

de s 



ds 



j'=i 



6 





dV A 






ds 2 


s=0 


and, because of symmetry, 




(m \A\ n) = — (m 


de sA 
ds 


s=0 



N 



n 



i=i 



(m \A ® A\ n) 



1 



N(N -I) 
Now we have the following result 



rn 



2„sA 



d e 



N 



ds 2 



s=0 



n 



m 



J2 A 



n 



(11) 
(12) 



(m | ®Are sA | n) 



-1/2 , N -1/2 



n 



max(n,m) 

x E 



TV 

x u x n x n x n 



((e sA hr ((e sA )nP l ((e sA )lt)^ T (V A )ll)* U , 



with m = x-ii + and n = x^ + x^ and A" = x^ + x-^ + x^ + x^. As seen in (11), 
we can calculate the derivate of the previous formula and divide by A" to obtain (n \ A\ m). 
This yields 

(m | A\ n) = (m 6 m>n A TT + y/m(N - m + l) <J m ,„_i A u 

+ y/(m- \){N -m)5 m ^ n A n + (N - m) 5 m<n A u ) . 



Similar computations with the second derivatives yield 

(m \B ®C\n) = m(m - 1) 5 TT C TT 6 m , n 

+ m y/m(N-m+l) (b tt C u + B n C n )<? ro _i, n 
+ m a/(A^ - m)(m + 1) ^ TT C iT + B U C^ 
+ m(N — m) i^B^Cii + B^C^ 



In particular, 



+ Vm(m - l)(jV - m + 2)(N - m + 1) B n C n 5 m - 2 ,n 
+ m(N-m) (B U C U + 5 iT C n )<5 

+ (N — m) ^m{N - m + 1) (s n C u + S u C n ) ^m— l,n 
+ \/(™ + 2)(m + l)(iV - m){N - m - 1) 5 iT C iT 5 m+2 ,„ 
+ (N — m — 1) y/(m+l)(N-m) (b u C u + 5 U Q T )<5 m+l,n 
+ (N-m)(N-m-l)B u C u 5 mjn ]. 

Ttn-2 I n)(n\ = ^(P 1 + P- 1 + P i + P_,) + O (1) 



where P e denotes the projection on (y/n I T) + e \/N — n | j) j for e = 1, —1, i or — i. 

We obtain that \n)(n\ is separable up to a correction of order jj. A similar computation 
shows that the marginal determined by (10) is, up to order l/N, separable. The following 
theorem provides a non-combinatorial answer to the question. 

Theorem 4. The distance between the two-site marginals of symmetric states on N sites 
and the exchangeable two-site states is not larger than d(d+l)/N where d is the dimension 
of the single-site algebra. 

Proof. Let us denote, for B G A^^(C), by Bj a copy of B at site j. By positivity and 
symmetry of an extension u N of uj we have 

N 2 

< cu n ((J2 B j) ) = N(N- l)u{B®B) +Nuu(B 2 ). 

Let P s be the projector on the symmetric subspace of C d <g> C d , which has dimension 
d(d+ l)/2, then 

Tr P s B <g> B = - Tr B 2 + - hi b) 2 for every B G A^(C). 

Choose now c = Nd(d + 1)/(N -l + d(d+ 1)) < d(d + 1) for d = 2, 3, . . . and N = 3, 4, . . ., 
then for every 5 G A4^(C) 



Tr P S P <g> 5 



> o. 

Now by theorem 2 we get that 

X G A* d (C) ® A^C) » (l - ^)c(X) + £ Tr P S X 

is an exchangeable state. And so we have that 

d(S~,S~)<^<*±±>. 



□ 



4 Mean-field models of composite particles 

The Hamiltonian of a mean-field system of iV quantum spins with a pair interaction h is 

HN = -^ E *Hr (13) 

{j,j'|l<i<j<iV} 



Here h is a Hermitian matrix on C d ® C d which is invariant under the flip operation 

(C ® 77 \h I y ® ^) = (v ® C |/i I ^ ® y?), j],(,^^eC rf . 

We shall, moreover, assume that h is ferromagnetic in the sense that there exist X a = 
(X a )* e M d (C) such that 

h = J2 xa ® xa - ( 14 ) 

a 

The factor 2/N in (13) is needed to obtain a good thermodynamic behaviour. 
A common example of such a model is the BCS-model where 

N . N N 

//»=-fc(i:$r)-^(E#)(Esr) 

«=i i=i j=i 

= - h & S ?)-2N ^ (S*S* + S?S«) +0(1). 
Here S^, S' 2 ' and S* 2 denote the generators of SU(2) 

and S ± = S x ±i S*. 

Using (14), we can rewrite the iV-particle Hamiltonian 

o j=l a j=l 

The second term in this expression has a norm of order 1 and is therefore thermodynam- 
ically irrelevant. Therefore, up to a correction of order 1, is a sum of negative terms. 
Moreover, the average ground state energy can essentially be computed by varying over the 
fully symmetric pure states, which is a proper subclass of the symmetric states, sometimes 
called the Bose symmetric states. 

As with exchangeable states, there is the notion of Bose exchangeable states. A state 
id on (gjjv-M^C) is called Bose exchangeable if it admits for any M > a Bose symmetric 
extension uj^m+m) to ®N+M-M-d(C). I.e., for any permutation rr of a set of N + M points 
and any A G <8>Ar+AfA4d(C) 

ld( N+M )(A) = id^ N+M ){AU^) (16) 

with Un as in (1). Note that the asymmetry in condition (16) is only apparent as 

id {N+M) (U n A) = id( N+M )(A*U n ) = id( N+M )(A*) = id {N+M) (A). 
The analogue of theorem 1 is then [7] 

Theorem 5. If id is a Bose exchangeable state on ®^M. d {C) , then 

id = d/i([y?]) <S>n M 
Jc d . 

proj 
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where Cp ro j is the complex projective d- dimensional Hilbert space and n is a probability 
measure on Cp ro j. By [ip] we denote the pure state of A4d(C) determined by the subspace 
Cip with \\ip\\ = 1, i.e. 

[<p](A):={<p\A\<p), AeM d (C). 

The asymptotic ground state energy density of a mean-field Hamiltonian with pair 
interaction h is then given by 

e (h) : = lim -j- iniu(H N ). 

AT— >oo iV cj 

Because of the permutation invariance and of condition (14), we have 

eo(/i) = -max ([<p]®[<p](h)). (17) 

Indeed, by theorem 1 it suffices to compute the infimum over product exchangeable states 
and if 

P = ^2n\<Pi){<Pi\ 

i 

is the eigenvalue decomposition of p we have, using condition (14) and the convexity of 
x I— > x 2 

p ® P (h) = e p( x °) 2 = E (E r * n (*«)) 2 

a a i 

<EE r «(NW) J = E r * N ® MCO- 

a i i 

The state space of a composite particle is of the form C dl <S> C d2 . We shall consider 
a simple pair interaction /i 12 = h± <g> h 2 between such pair interactions with ft^ and /i 2 
ferromagnetic in the sense of (14). We now have the following result 

Theorem 6. Assume that hi e .M^C) ® .M^C), i—1,2 are Hermitian, invariant under 
the flip and satisfy condition (14)- Assume, moreover, that hi is positive definite, then 

e (/ii <E> h 2 ) = -e (/ii) e (/i 2 ). 

Proof. By the negativity of the mean-field Hamiltonians corresponding to pair-interactions 
satisfying (14), see (15), we have 



e (hi <S> h 2 ) — - max ( [ip 12 \ <S> [<pi2](hi <S> h 2 ) ) 

[V12] V / 

< - max ( [(p 12 \ ® [ifn] (h <g> h 2 ) ) 

{[V12]|[¥'12] = [¥'l]®[¥'2]} V / 

= - max ® [</?i](/ii)j max ([</? 2 ] <g) [</? 2 ](/i 2 )) 



= -e (hi) e (h 2 ). 

To obtain the converse inequality, consider a normalized vector ip 12 e C dl <8> C ' 2 and the 
state 

U2 W®bi 2 ](/n®i) 
10 



on M. d2 {£) ® M. d2 (C). This state is flip-invariant and, because 



h 1 = ^X a ®X a 

a 

enjoys the property 

J? 12] {Y ®Y) > 0, Y = Y* e M d2 (C). 
Hence, by theorem 2, it is a mixture of product states. Then by the remarks above 

ui [ r 2 \h 2 )>e (h 2 ). 

We therefore have 

e (/ii ® h 2 ) — - max ( [ip 12 ] <g) [y?i 2 ] (hi ® h 2 ) ) 

= - max ([ip 12 ] <g> [^i 2 ](/ii <g> 1) ^ 12l (/i 2 )) 
[V12] v / 

> e (/i 2 ) max ( [(p 12 \ <g) [</?i 2 ](/ii (8) 1) ) 

> -e (/ii) e (/i 2 )- 

The last estimate follows from the fact that 1 is positive definite and satisfies condition (14). 

□ 

Two remarks are here in order. There doesn't seem to be a simple extension of theorem 6 
to finite temperatures, at least no simple relation between the free energy densities of the 
composite system and the components seems to exist. A second remark is that the theorem 
can be used to give a partial answer to the problem of multiplicativity of maximal 2-norm 
of quantum channels [4, 8]. Unfortunately, the positivity condition on hi imposes some 
restriction on the allowed channels. A further elaboration of this matter will be considered 
in a future publication. 

Appendix A: The map Vd 

Every Hermitian matrix B in Ai^(C) can be written as 

R _fb (V 

where b e R, | ip) is a vector in C d_1 and B a matrix in A / fJj_ 1 (C). We then define the map 
V d : M d {C) -> R d2 inductively as 



V d {B) := 



/ b \ 

v^Rel^) 
\/2Imj ip) 
\V d -i(B )j 



This map has the following properties for B 1: B 2 e M d (C) 
i) V d (B 1 + B 2 ) = V d (B 1 ) + V d (B 2 ). 
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ii) For every A G R, V d (XB 1 ) = AV^Bi). 

in) TrB 1 B 2 = {V d (B 1 )\V d (B 2 )). 

This can easily be proved by induction on d. Moreover, the map V d is one-to-one and onto. 
Note, however, that the map V d is basis dependent. 

Appendix B: The map 

The subspace 7C 

Before we start to search for a good map M d , we take a closer look at the subset /C 
of flip- symmetric, complex, hermitian matrices on C d2 . We begin by decomposing the d- 
dimensional Hilbert space C d in a direct sum of a one-dimensional and a (d— l)-dimensional 
space, C d = C © C d_1 . We are interested in the symmetric, (C d <g) C d ) s , and antisymmetric, 
(C d <S> C d ) a , subspaces of C d <g> C d as they are the ones left invariant by the elements in C*. 
We consider a basis {e , . . . , e d -i} of C d . Then a basis of (C d <S> C d ) s is given by 

{e <g> e , #1, ■ ■ ■ , gd-i, fi, ■ ■ ■ , fd(d-i)/2} 

where (?j := ^(^o <E> + e« ® eo) and where the /j generate the symmetric subspace of 
C d_1 <g) C"' -1 . Similarly, a basis of (C d ® C d ) a is given by 

{hi, . . . , h d _ u hi, . . . , fc(d-2)(d-i)/2} 
where /ij := ^75 (eo <8> — e« <8> eo) and where the ki generate the antisymmetric subspace of 

A matrix A G /C can be written in this symmetric-antisymmetric basis as 

/a (y> | ($ | \ 
(p) X x Fi 

$) Zi o o (is) 

X 2 Y 2 

\o o or; z 2 ) 

where 

a G C, e C d ~\ $ G (C^ 1 <g> C^ 1 ) 8 , X 1 , X 2 G Mi_i(C) 

Fi : (C d_1 (8) C^ 1 ) 8 -> C d_1 and Y 2 : (C d_1 ® C d_1 ) a -> C d_1 . 

In order to ensure that we map the subspace /C in a suitable vector space, we can 
count its real dimension. The restriction of elements of /C to the symmetric subspace needs 
d(d + l)/2 real parameters on the diagonal and 2 times (for the real and imaginary parts) 
[d(d + l)/2] l(d(d + l)/2) — l]/2 off the diagonal. For the restriction to the antisymmetric 
subspace we need d(d — l)/2) + [d(d — l)/2][(d(d — l)/2) — 1] parameters. In total this 
amounts to d 2 [d 2 + l)/2 real parameters, which is exactly equal to the dimension of the 
symmetric real matrices of dimension d 2 , i.e. the matrices M G Ai d 2(R) such that M = M T 
where T denotes transposition. 
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The map M d 

Denote the symmetric real matrices of dimension d 2 by .Mjj 2 ([R). Using the parametrisa- 
tion (18) for A E K, we define the map M d : IC ^ M#(R) by 



M d (A) := 



a 

| Rey?) 
(X X + X 2 ^ 



(Re V | 
ReXi - ReX 2 



[Re*] 



(limp | 
ImXi — ImX 2 



-X 1 +X 2 - 



[Im*] 



'ImXi-ImXa |FT A1 \T ReXi - RcX 2 
+ [Im*] [Re*] 



Ti(Yi,y 2 ) 



Ti(yi,y 2 ) T 



T 2 (Yi,y 2 ) T 



M, 



T 2 {Y U Y 2 ) 
'Zi 

z 2 ;y 



where for 2 7^ j 

[Re$]ji := Re^, [Re<% 



-j=Ete$ij, [Im$] u := Im<5>ii and [Im<% := -^Im^. 



We describe the maps Ti and T 2 in the two following paragraphs. As with V d , the map M d 
is basis dependent 

The map Ti Recalling that {ej}^"/ is a basis we choose in C d_1 , let us, for i < j, 
i,j — l,...,d—l and any matrix B Q G yWjj_ 1 (C) put 

f3 R (i,j) := a if and only if (V d -i(B ) \e a ) = V2Be[B ]ij 

:= a if and only if {V d _i{B Q ) \e a ) = \Z2Im[B ] i:j and 
:= a if and only if (V d -i(B ) \e a ) = [B ] u . 

This way of denoting the matrix elements will be useful later on when we will compare 
B <S> B with the projection on V d (B). We also define e e k = 1 if k < £ and —1 otherwise. 
We are now ready to define the map Ti by looking at each of the matrix elements. In the 
following, i, k, £ run from 1 to d — 1 and i<£, i^k, £^k 



• [Ti(Yi,r 2 ) 

• [Ti(Yi,r 2 ) 

• [Ti(Yi,r 2 ) 

• [Ti(Y u Y 2 ) 

• [Ti(Yi,y 2 ) 

• [Ti{Y u Y 2 ) 

• [Ti{Y^Y 2 ) 

• [Ti(y!,y 2 ) 



i,/3(i) '■— Re[^l]i,ii 

fc,/3(i) := -j2=Re([Yi\i,ik + 4 P^kifc) 



^Re([F 1 ]^ + 



+ 



[Y^u+_4 [ Y 2]e,u \ 
V2 ) 



k,(3 R (i,. 



e,/3i(i,e) '■- 



k,/3i(i,£) 



Re 



\Yi]e,ik + 4 [Y2\e,ik + \Yi]i& + 4 



[YiU + 4 m 



^Im([r 1 ]^-^ + e ^ F2k 



-Im 



V2 
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The map T 2 The notations are similar to the ones used for the map 7\. Again we define 
each matrix element 



• [T 2 (Yi,r 2 ) 

• [T2(Y U Y 2 ) 

• [T 2 (Y 1 ,Y 2 ) 

• [T 2 (F!,F 2 ) 

• [T 2 (Y U Y 2 ) 

• [T 2 (Y 1 ,Y 2 ) 

• \T2{Y U Y 2 ) 

• [T 2 (F!,F 2 ) 



i,/9(i) : = I m [^l]i,M 

k,p(i) ■= ^Im^[Yi] i>ifc + e l fc [y 2 ]j,ife 

[nk + 4 g^Mg 



i,p R (i,t) ■= 75 Im^i]^ + ^ 



i,Pz(i,e) '■- 



1 nifyl [frl]^ + 4 [^j^A 

^Re^Yij^ -j= j 

-73 Re (J r ik« 7f J 

/ ftlkifc + 4 P^k.tfc - [Yl]j/k ~ 4 [ y 2]»,£fc 

V 2 



:= Re 



)■ 



One can easily see that, given Ti(Yi,Y 2 ) and T 2 (Yi,Y 2 ), one can reconstruct the matrices 
Yi and Y 2 . Also these two maps are real linear. 

Properties of the map M d 

The map Md has similar properties as the map Vd 

• M d {A l + A 2 ) = M d {A 1 ) + M d (A 2 ). 

• For every A G R, M d (AA) = AM d (A). 

It is also one-to-one and onto. Again one can easily check these properties by induction on 
d using lm(X[ j - X$) = -Im(Xf - Xf ). 



The image of B ® B 

Fix i? G .M^(C) and consider the tensor product of B with itself 

/ b 2 V2b(ip\ ((^®^) s | 



B®B = 



\ 





y/2b\il>) ftSo+lVXV'l ((^l®5 ) s 

|(^®^) s ) (|^)®^o) s (B ®B ) S 

fcSo-lVXV'l «^|®5o) 

V (|V)®5 ) a {B ®B )*J 



We will prove that 5 <g> B is mapped by M d on 1 1^(5)) (1^(5) | with 

/ b \ 

V2(Reip | 



V d (B) 
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First we write down the image of B <g> B 




bV2\Reip) Re| ip){ip \ + [Rc(V> ® V) s ] Im| V} (V> I + [M^ ® VOI 

6v / 2|In«/'> Im| V}(V> I + [Im(V"X> V>) s ]* Re| V>W> I - [Re(i/> ® </>) s ] 

VHV d -i(Bo)> T 1 (((i/.|®B ) s ,((V|®B ) a )* T 2 («V|® Bo) s ,((V>l®Bo) a )* 



&<V d -i(Bo)| 
Ti(((V|® Bo) s ,((^|® B ) a ) 

M d _!(Bo®Bo) 



The first row and column are encouraging but we still have some steps to verify. If we 
use induction on d, we also get that M d _ 1 (B <g> B ) = | Vd-i(B )) (V d _ 1 (B ) |. Let's look at 
the other parts of the matrix. 

Looking at the elements in the middle of the matrices M^iB <g> B), we need to prove 
that 

• Re| VX^I + [ReO®^) s ] = | v / 2Re^)(v / 2Re^ |, 

• Re | i/>){i/>\- [Re{ip <g> ip) s ] = | v / 2Im^)( v / 2Im^ | and 

• Im| ij)) (ijj | + [Im(V> <g> ^) s ] = I Re^) (v 7 ^ Iim/> | 

in order to obtain that B <g> B is mapped on | V(B))(V(B) \ . 
Let's look at the different matrix elements 

• Re | ip){ip\ + \Re(ip <g> vp) s ] = 2| Re'0)(Re'0 |. Indeed, it is easy to see that 

[Re| I + [Re(V> <g> V) s ]]u = ((Re^) 2 + (Im^) 2 ) + Re^ 2 = 2(Re^) 2 and 

[Re| i/>)(i/>\ + [Re(V> <S> ^) s ]]ij = Re^i Reijjj + Im^ Im^ + — Re(v / 2^ 

v2 

= 2 ReV'i ReV'j • 

• Re | i>){i> | — [Re(^ <g> ?/>) s ] = 2 | Inr0)(lnr0 |. The proof is similar to the one above. 

• Im| ip)(ip | + [Im(^ <S> ip) s ] = | Re*0)(lm0 |. Indeed, 

[Im| + [Im(^ <S> ?p) s }}ii = 2 Re^jlm^j and 

[Im| ip){ip | + [Im(^ ® V ; ) s ]]ij — ReV'j Im^- — Im^j Re^- + -^= lm(V2ipi ipj) 

V2 

= 2 Re0j Im^j . 

Now the proof is almost complete. We still have to verify that that ((01® B ) s and 
{{ip | <g> A,) 51 are mapped by Ti and T 2 on | v^Re^^-i^o) | and | V% Imip) (Vd_i(-B ) | 
respectively. 

The map Ti We now verify that 



• | ® B ) B ,((^ | ® £o) a )]i,/j(i) = Re(v / 2^ [B ]«) = Re& [B ] 

•[^i((^l®^o) B ,(^|®5o) a )] fcl/3 ( 



15 



= ^Re(V>i [B ] ik + ij k [B ]ii + i k e\ [B ] ik - ip k [£„]«)) 
= V2Re*p k [B ] u 

.[TM^\®B y^\®B r)] hMt/) 

1 r> i AT / r d i , ^ t^oltf + ^ [B ]u + (A [Bq] u - ip e [_B ]u) 
= -^He{V2ipi [B \ ei H -j= ) 

= Re(^([So]« + [£,,]«)) = V2Re^ V2 Re [5 k 

1 D / AT , r n 1 , ^ t^oltt + ^ [-B ]« - (^i [B ]« - ipe [B ]ei) , 
= ^*te\y*V>t [Bq\u H ^= j 

= Re(M[BoU + [Bo]u)) = V2 Re^ V2 Re[B ] it 

= ^Re(tpi [B ] ek + i) k [B \ u + e\ e\(^i [B ] £k - ip k [B ] ei ) 

+ if> t [B ] ik + i) k [B ]u + e[ e k e (ip t [B ] ik - tp k [B ]u)) 
= Re^ fc ([S ]tf + [B ]«) = v^Re^ >/2Re[B ]tf 

= -^ im l-V2^ [i^ojtt H ^= j 

= Im(^([S ]tf - [S ]«)) = v^Re^ v^Im[S ] tf 

= — Im(v^^ [fl ]tf - — ^ U + — ^* ~ — ~ — ) 

= Im(^^([5 ] !f - [£„]«)) = Re^ >/2Im[£ ]tf 
•[ri(«^l®5o) B ,(^|®5 ) a )] fcl/3j ( il0 

= 2 Im ( - ^ - ^fc [-B ]« - ef(^ [-B U - i) k [B ] ti ) 

+ ip t [B ] ik + i) k [B ]u + e e k [B ] ik - i\) k [B ] u )) 
= Im^fc([S ]tf - [S ]«) = V2Re*p k V21m[B ] u 

The map T 2 The proof that 

T 2 (((^ | <g> Bof, | ® S ) a )) = | v / 2Im^)(Vd_ 1 ( J B ) | 

is completely similar, so we will not provide the details. We have now proven a one-to- 
one correspondence between B <g> B e (•Md(C) <E> -Md(C)) h and the subset of rank one 
projections in .A/f^lR). We now have real linear one-to-one and onto maps V d and M d that 
satisfy condition (Hi) Section 2.2. Let us now examine condition (ii). 

Appendix C: Tr A (B®B) = (V d (B) \M d (A)\V d (B)) 

We start by calculating the trace of A(B ® B). 
TiA(B®B) = TiA s (B® B) s + Tr A a (B® B) a 
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a b 2 + bTr B X 1 + (ip \Xi\ ij)) + Tr Z x (B <g> B f + 2 ReV2 b (if> \<p) 
+ 2 Re(^ <g> i) |$) + 2 ReTr((^ | <g> 5) s Y^J + \b Tr fi X 2 - (V> |X 2 | V) 
+ Tr(B ® 5 ) a Z 2 + 2 Re Tr((0 | <g> £ ) a Y~ 2 *] . 



We can restructure this expression 



Tr A(B ® B) =bab + 2 bV2Re((p \ip) + 2b Tr( Xl + X<1 ) B 

+ {ij; \X, -X 2 \iP) + 2 Re(V> ® V 1$) 
+ 2 Re Tr((^ | <g> Bf Y* + 2 Re Tr((0 | <g> 5 ) a F 2 * 
% 



+ Tr(£ <g> B ) 



Z, 



We rewrite the first line of the right-hand side of the above equality. To make the link 
with Vd(B) and M d (A), we express ip and ip in their real and imaginary parts. We also use 
property (in) of the map Vd- We then get 



bab + 2bV2Re(<p\iP) + 26Tr( 1 2 ) B 

= 5a6 + 26((Re V 9|v / 2Re^) + (Imp) [v^W) + 2 b (v d ^ ( Xl ± | V d ^(B Q ) 

This looks promising, we can also try to express the second line in term of elements 
appearing in Vd(B) and M d (A) or by looking at the real and imaginary part of the matrix- 
and vector components 

(ij \X X - X 2 \ ij) + 2 Re((V> ® |$) 

= ^((Re^) 2 + (Im^) 2 ) {[X&i - [X 2 ] u ) 



+ 2 [ ( Re ^ Re ^' + Im ^ Im^) Re ( [X^ - [X 2 ]ij) 

{hj\i<j} 

- (Retpilmtpj + Im'0 i Re^j)lni([Xi] ij - [X 2 ]y)j 
+ 2 [((Re^) 2 - (Im^) 2 )Re$ii + 2 Re^i Im$. 



+ 2 ^ (Re^i ReV'j - Im^ Im^Re^- + (Reipi haipj + Im0i Re^)lm<l> 

{»,j'I»<j} 

ReX : - ReX 2 



= ^Re-0 
+ ^\/2Im^ 
+ 2/v / 2Re^ 



+ [Rc$] v^Re^ 



ReXi - ReX 2 



ImXi — ImX 2 



- [Re$] V2hmp 
+ [Im$] v^Im^ 



This also points out to the equality we are trying to prove. The fourth line is less straight- 
forward but we can rewrite it 



2 Re Tr((V> | <8> B) s Y* + 2 Re Tr((0 | <g> £ ) a Y* 
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= 2Re [J2 ^ [ F ik(fc,fc)^fc [B ] ik 

i k 

+ { ^ B ^ i£ + ^ l i B o\ik) + tyk [B ]u - [B ]ik) \^2\i,{k,t) }] 

{k,e\k<£) 

= 2 Re [v^ [Yi}i,(u)4>i [Boh + Y ^ Fik(M)^ [B }ik 

i {k\k^i} 

+ y M^u(\yiWi,k) + 4^2]^)) 
{k\k^} 

+ Y ^ [ s o]tf([^i]i,(/,fc) + 4[^}i,(e,k)) 
{k,e\k^e, e^i} 

= 2Re^[v / 2[F 1 ] a ^+ ^ ^fc(Pi]i 1 «,fc) + 4[^2]i ) (i 1 fc))[So]i i 

+ ^ (v / 2[Fik( V) ?/Y+ ^ ^([^i]i,(/ 1 fc) + 4[^]i ) (/ 1 fc)))[So]« 
{WO 

= 2Re^[v / 2[Fi] a ^+ ^ ^fc(Pi]i 1 «,fc) + 4[^2]i ) (i 1 fc))[So]i i 

i WMO 

+ (v / 2[?i],,(V)^+ ^ ^([^ik(^) + 4[^]i 1 (< ) fc)))[S ]« 

{i\i<£} {k\k^£} 

+ (v / 2[n]^,M^+ Y ^([ ? ik«,fc) + 4[^k(i ) fc)))[So]tf 
{^|j<^} fflfc^O 

= 2 ^[(v^Refy^,^) Re^ + v^Im^]^ Im^) 

+ ( Re ^ fc Re ([ y i]i,(»,fc) + 4 P^li.Ci,*)) + Im^fc Im([Yi] ii( j >fc ) 

{fc|fc#0 

+ 4 [*2k(i,fc))) [-Bolii 

+ S { (v^Re^ili.^Re^ + v^ImfYil^^Im^ + Rep^^Re^ 

{i\i<e} 

+ v / 2Im[y 1 ] A(i)i) ImV'^ 

+ J] ( Re ^ fc Re ([ F ik(^) + 4 [*2]i,(<,fc)) + Wfc Im([Yi]j ; (^ fc ) + e£ [^k^fc))) 
{fc|fe^} 

+ X (^Re^fc Re([Fi] £i (i ifc ) + 4 [^k^fc)) + Im0i Im([Yi] £>( ; ifc ) 
{fc|fe^i} 

+ 4 [^]£,(, fc )))Re[ J B ]^} 
+ X { ("V^ Re[Fi] av) Im^ + v^MY^^Re^ + v 7 ^ Re[Fi],, M Im^ 
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+ ("Wfc Re([Yi]^ )jfc) + e e k [Y 2 } um ) + Re^ fc Im([Yl]^ ifc) + e e k [Y 2 ] um )^ 



+ ^ Re([^i]<,(i,fc) + 4 P^kftk)) - Re^fc I m ([^i]^,(z,fe) 

{k\k^i} 

+ el[Y 2 ] Uhk) ))lm[B 



= 2^v / 2Re^ T X {J X ,Y 2 ) ^(Sq)) + 2 (v^W T^Y,) V^_i(S ) 
Finally, using induction on d 

% 



Tr(5o ® .Bq) 



^-l(S ) 



Zi 

z 2 



z 2 

If we put this all together we get what we wanted to prove, namely 

TrA(B®B) 

= 6a6 + 26((Rev?|v / 2Re^) + (limp | V2 Imip)) 
'X 1 + X 2 



^_i(So) 



+ 26^1^ 



V^-i(So) 



+ ^Re^ 
+ 2/ v / 2Re^ 



ReX : - ReX 2 



ReXi - ReX 2 



+ [Re$] v^Re^ 



- [Re$] v^Im^ 



ImXi — ImX 2 



+ [Im$] v^Im^ 



+ 2 (v^Re^ + y 2 )| V d ^(B )) + 2 (v^W (T^ + Y 2 )\ V d ^(B )) 



+ (V d ^(B )\M d 



Z 1 



z 2 

(V d (B)\M d (A)\V d (B)). 



)\V^i(Bo)) 



To summarize, we have found maps 

V d : M h d (Q -> R d2 and M d : (.M^C) <g> A^ d (C)) h -> A^(R) 
with properties that allow us to prove the second part of theorem 2, see section 2.2. 
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